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Abstract 

Recent studies suggest that PML (Promyelocytic Leukemia) suppresses p53 protein degradation by in¬ 
hibiting MDM2 protein in the nucleus [T] , and regulates a number of biological functions [5]. We modeled 
a PML-MDM2-p53 regulatory network by integrating p53-MDM2 autoregulatory model [3] and related 
PML pathways. We found that p53 dynamics switched at various dynamical states induced by PML 
which corresponds to different cellular states. Our results show clear transitions among the stabilized, 
damped and sustain states under different stress conditions induced by PML in the system. These states 
in p53 dynamics is the signature of existence of various cellular states and the phase transition like 
behaviour of these states involve various biological significance. 

Keywords: PML (Promyelocytic Leukemia), PML-MDM2-p53 network. Phase transition. Determin¬ 
istic model. 


Introduction 

PML is a well known tumor suppressor nuclear protein [I]. It has a small nuclear structure and present 
in almost all mammalian cells. It regulates many processes such as cell division, terminal differentiation 
of myeloid precursosr cells and neural progenitor cell differentiation [4]. It also acts as a transcription 
factor for many important regulatory proteins which helps in initiation of DNA damage response, DNA 
repair mechanism and organization of chromatin. A number of regulatory proteins are associated with 
PML such as SUMOl (Small ubiquitin-related modifier 1), SUM02 (Small ubiquitin-related modifier 2), 
MDM2(Mouse double minute 2), Tp53(tumor protein p53), CREBBP (CREB binding protein), RARA 
(retinoic acid receptor, alpha), ZFYVE9(zinc finger, FYVE domain containing 9), UBC(ubiquitin C), 
DAXX(death-domain associated protein), SPIOO (SPIOO nuclear antigen) etc [IHZ]- PML is a key protein 
for generation of PML-nuclear bodies (which is also known as kremer bodies), nuclear domain 10 or PML 
oncogenic domains [5]. PML is very sensitive as well as responsive to various stress like DNA damage, 
viral infection, transformation, and oxidative stress etc [5]. Recent studies suggest that PML can suppress 
p53 tumour suppressor protein degradation by controlling MDM2 protein in the nucleus [I] . 

p53 is a widely studied and an important protein in study of normal as well cancerous cells. It 
controls many important regulatory pathways such as cell cycle, apoptosis, DNA repair, cell differentiation 
etc [aiiniin]. It is responsive to the various stress like osmotic stress, oxidative stress, DNA damage 
by UV and IR, chromosomal abbration, chemical exposure, hormonal variability etc. In recent research 
reports various models have been proposed to shown the p53 activation through different stresses m- 
MDM2 is a oncogenic protein which acts as negative regulator of p53 protein. MDM2 protein acts as 
E3 ubiquitin ligase and voluntarily interact with p53 and forms p53-MDM2 complex and this leads to 
decrease in p53 concentration level within the cell [HE]. Several studies have been performed to show 
the p53 and MDM2 interaction and control mechanism of p53-MDM2 network. Sevaral attempt are done 
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to understand the mechanism to control p53 inside the cell but still there is scope to understand the p53 
restoration and the molecular event associted within the system. 

In the present work we designed and studied a p53-MDM2-PML integrated network by incorporating 
PML to the p53-Mdm2 regulatory network. In this model PML protein interacts with both p53 as well 
MDM2 protein and activates p53 and induces p53 expression within the cell which leads to increase in 
p53 concentration level within the cell. Several experimental and theoretical models have been proposed 
to study the PML interaction with p53-MDM2 network but still the temporal bahaviour and dynamical 
states of p53 induced by PML is not fully understood. The role of noise in the system is not systematically 
incorporated and investigated. We study some of these problems in this work in order to understand 
how does PML interferes cellular activities via p53. Our work is organized as follows, PML induced p53- 
MDM2 biochemical regulatory model and its molecular interaction is described in section 2, numerical 
results are presented with discussions in section 3 and some conclusions are drawn based on the results 
we obtained in section 4. 


Materials and Methods 

PML induced p53 — Mdm2 regulatory model 

p53-MDM2 autoregulatory process is considered to be a very important regulatory process since it in¬ 
terferes most of the cellular activities both in normal as well as cancer cells mM- p53 is reported to be 
synthesized in cell with a rate constant [5]. Mdm2 is a negative regulator of p53, and its translation 
from Mdm2_mRNA is considered to be occured with a rate ki. p53 acts as transcription factor which 
interacts with Mdm2 gene due to which transcription of Mdm2 gene into Mdm2jrnRNA occurs with a 
rate ^ 2 - The half life of Mdm2jmRNA is low which leads to the degradation of the Mdm2jmRNA which 
is considered to be occured with a rate constant fea. The ubiqutination of the Mdm2 protein is taken to 
occur with a rate constant k^. Mdm2 interacts with p53 protein in nucleus and forms pbZ-Mdm2 complex 
with a rate constant fcy. Further, Mdm2 which act as ubiqiutin ligase which leads to the degradation of 
p53 concentration in cell with a rate fcg [11]. The dissociation of the pbiJkldm2 complex is considered 
to be occured with a rate fcg. PML is considered as a nucleoprotein which is constantly expressed within 
the nucleus and varies in system to system with rate kpML- PML is reported to interact with Mdm2 to 
form PML_Mdm2 complex with a rate constant feio. Furhter the degradation of Mdm2 occurs due to 
its interaction with PML, and is assumed to occur with a rate constant ku. Moreover, PML changes the 
structural configuration of Mdm2 so that it is not available to interact with p53 and due to which p53 
concentration increases in the cell. PML also directly interacts with p53 and forms PMLjp53 complex 
with a rate constant ki^. The interaction of PML and p53 leads to the phosphorylation of p53 protein 
on serine 20 [HHn]. This phosphorylation of p53 saves it from ubiquitation induced by Mdm2. Further, 
the dissociation of PMLjpbS complex takes place with a rate constant ku. Moreover, degradation of 
the PML protein due to its shorter half life is reported to be taken place with a rate consrtant ki 2 |9|. 
In Table 1, we listed molecular species associated with p53-MDM2-PML network. In Table 2, we listed 
the reaction channels, their rate constant and values taken in the simulation. The biochemical reaction 
network of the model is shown in Fig. I. 

Let the configurational state vector of the system we consider corresponding to state variables (molec¬ 
ular species in Table 1) be given by x(t) = {xi,X 2 , ...,a;„}^, where, T is the transpose of the vector and 
n = 7 Based on the biochemical network shown in figure I, we have translated the chemical reactions 
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(Table 2) into a set of coupled ordinary differential equation using mass action law given by, 


dxi 

dt 
dX2 
dt 
dX3 
dt 
dx4 
dt 
da; 5 
dt 

dxQ 

dt 
da; 7 
dt 


- k^XlX2 + k^x^ - ki^xiXr, + kl4X^ ( 1 ) 

fcia;4 - k4X2 + kQX3 - krXiX2 + k^x^ - kiQX3X2 (2) 

—kexs + k7XiX2 — ksxs ( 3 ) 

k2Xi — k^Xi (4) 

kpML ~ kioXc,X2 + kiiXQ — ki2X3 — ki^XiX^ 

+ki4Xr - kiQx^xe (5) 

fcioa;5a;2 - fciia^e (6) 

fci3a;5a;i - kiAX7 (7) 


where, {ki} i = 1,2,... ,M{M = 14) represent the sets of rate constants of the reactions listed in Table 
2 and names of the molecular species listed in Table 1. We have used runge kutta 4th order algorithm 
for numerically solving the set of differential equations of our system. We developed our own coding 
simulation software which is written in fortran language [18] . 


Results and Discussion 

We present our numerical results as well as discussion on the results. We also tried to relate our results 
with possible biological processes arised due to PML interaction with p53-Mdm2 regulatory network. 

p53 — Mdm2 bio-chemical network activation 

The temporal behaviour of p53 protein for different values of PML creation rate constant kpML (Fig. 
2) shows different states driven by PML. The small values of kpML ((0.0008) could not able to activate 
from its stabilized state which corresponds to nearly normal state (Fig. 2 lowermost panel). However, 
for slight larger values of fcpMi ((0.002), p53 dynamics exhibits damped oscillation for certain range 
of time (activated state) which depends on kpML and then maintains stabilized (come back to normal 
condition). Further, increase in kpML elongates the time of activation, and p53 dynamics become sustain 
oscillation state (~ 0.006) showing strongest activation induced in the system. Further increased in kpML 
forces the p53 dynamics to damped oscillation and then reached stabilized state at larger concentration 
level (signature of apoptosis). Further, it is observed as the concentration level of PML increases which 
corresponds to the rate constant of synthesis of PML within the cell, the temporal behaviour of p53 is also 
changes. This study suggests that as PML concentration increases the stress level in the cell also increases 
which is suppose to be due to molecular interaction of PML with p53-Mdm2 network. Further, It is also 
to be noted that as value of fcpMi increases from this particular range, the amplitude is also increases. 
This is probably due to higher stress generated by PML, which induces the p53 activation. This indicates 
that at this particular concentration of PML, the system is succumb to the stress and system is moving 
towards the apoptosis(cell death). Finally, when the concentration of PML is reached to be very high i.e 
at fcpML=0.05, we have got steady state level of p53. This suggest that, the system is succumb to stress 
and there is no further activation of p53 is observed and and cell goes for apoptosis(apoptosis) [T9] . 

Similarly, we have analysed the temporal behaviour of Mdm2 protein as shown in Fig. 3 for same 
values of kpML as in the case of p53. Similar effect is observed but correspodingly in reverse manner as 


4 


shown in p53 temporal behaviour which is supposed to be due to the sequestering effect of PML upon 
Mdm2 which is shown in Fig. 3. 

A two dimensional plot of p53-Mdm2 network and its associated molecule is also shown in Fig. 4. In 
the first left hand Panel in Fig 3, the quantitave measurement of p53 verses Mdm2 which shown through 
correspoding configuration space for various values of kpML- These results further verify the temporal 
behaviour discussed above. 

Stability analysis of p53 and Mdm2 

Stability analysis for the system is shown in Fig. 5 and 6. We have analysed the stability of the system by 
two parameter, namely amplitude and maxima in the dynamics (p53 and Mdm2). In Fig. 5 the maxima 
of p53 and Mdm2 are found to be remain stable up to certain values of fcpMi, then it increases with 
respect to kpML values up to 0.006 and after which it decreases i.e in between fcpML=0.006 to 0.008. 
Further increasing fcpMi values become constant. This observation suggests that as PML concentration 
increases within the systems, the stress upon the system increases which leads to the increase in synthesis 
of p53 in the system which we can easily understand from the increase in maxima in p53 and Mdm2 and 
corresponding amplitudes (Fig. 6). We can also observe the amplitude death regimes (for small and large 
kpML values), and in between stress or activated regime. 

Steady State Solutions for the systems 

The steady state solutions of the deterministic equations given by (l)-(7) corresponding to the concen¬ 
trations of the molecular species can be obtained by taking, d[xi]/dt = 0, where i = 1,2,..., N {N = 7). 
Imposing these conditions we get seven differential equations involving seven variables corresponding to 
concentrations of the seven molecular species. Using these equations the steady state solution for bi], 
[^lldeter can be obtained and is given below. 


i^^ldeter 


I k ^kjk^ 

ik2kr 


h 

ks 


( 8 ) 


This steady state solution gives us the behaviour of the p53 protein concentration when it is in 

extreme conditions and it is affected by various mechanisms. The equation (|5]) indicates that p53 level is 
in fact influenced by various rate constants of various molecular species, their complex formation, decay 
and creation etc. If we make the rate of formation of MDM2, ki a variable then from this equation it is 
found that is inversely proportional to indicating the inhibiting activity of MDM2 to the 

steady state level of p53. The factor fcs/fce can be seen as the rate of dissociation of p53_MDM2 per p53 
decay rate i.e. a factor which indicate average amount of availability of p53 and MDM2 in the system. 
As fcg/fce increases more p53 and MDM2 are available in the system and vice versa. 

We then solve for steady state solution of MDM2 i.e. [x2]*deter explained above and is given by. 


\X2\(leter 


Ikik2h A ^ 

k^k^ikj V kQ 


(9) 


Similarly, [x2]*deter gives us the steady state level of MDM2 in extreme conditions and as in the case of 
p53, its steady state level is affected by other rate constants. However this steady state level is directly 
proportional to the which is the rate of p53 protein formation. This means that the MDM2 steady 
state is strongly activated by p53 protein interaction in the network which is in agreement with the 
experimental results reported in qualitative sense. 
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Similarly, following the same procedure the steady state solutions of the variables and [ 2 ^ 4 ] 

are obtained by solving the steady state equations and are given by, 


[ 2 ^ 3 ] 


* _ 

deter 


k2k4:h A h\ 
kik^kr \ keJ 


( 10 ) 


NSeter = r 

From equation (fTOl) that [2;3]2gter i'^dicating suppressing effect of MDM2 on MDM2_mRNA, 

however, [a^aj^g^g^ oc which reveals that p53 activates MDM2_mRNA. Further, from equation (ITTl) 
we have oc fcs indicating p53_MDM2 is activated by p53. Now we solve for by relating to 

the rate constants of p53 and MDM2 as in the following, 


N 

[xe] 


deter 


deter 


fell k^kikr 

(1 + U) 


( 12 ) 


At constant level of [a^ej^eter; ^^^e inversely proportional of [a^sj^gtgg to gives rise the suppressing 
activity of MDM2 to PML steady state level which is in agreement with the experimental results and 
predictions. 


Stress driven by PML and Mdm2 interaction 

PML interacts with p53 and Mdm2 allowing to change the state of the system. The interaction of PML 
and p53 is one important interaction which can change the fate of the system via p53. The interaction rate 
of PML and Mdm2, fcio indicates the change in Mdm2 due to change in PML concentration level. The 
increase in fcio allows the increase in amplitude of p53 as well as Mdm2 indicating increase in activation 
in the system (Fig. 7). Further, it is also observe that as kio increases the amplitude death regime in 
small kpML regime (normal state regime) become shifted towards smaller kpML direction. It means that 
increase in kiQ the system is activated quicker. Similarly, the increase in fcio extends the amplitude death 
regime. Surprisingly, this increase in kiQ needs larger values of kpML to reach apoptosis. Same scenario 
is obtained for Mdm2 case also. 


Conclusion 

PML protein is found to be a tumour suppressor protein and it activates many protein to acheive cellular 
apoptosis. Although it directly interacts with both p53 and Mdm2, its behaviour is different for both the 
proteins. In the first case it activates the expression and in second case it modified the protein properties 
so that it can not negatively feedback upon previous one as reported in various experimental results. 
Our results support the stabilization effect of PML upon p53 and Mdm2. Simulation results clearly 
explain the temporal behaviour of the system with respect to stress driving the system to different states, 
stabilized, damped and sustain oscillations. 

The interaction of PML and Mdm2 enhances the activation of the system and resists the system to go 
to apoptosis. Similarly, interaction of PML with other interaction partners may influence in the cellular 
stress. The increase in PML level in cellular system may trigger cancer in some cases which is needed to 
be investigated further. 
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Table 1 - List of molecular species 


S.No. 

Species Name 

Description 

Notation 

1. 

p53 

Unbounded p53 protein 

Xi 

2. 

Mdm2 

Unbounded Mdm2 protein 

X2 

3. 

Mdm2 p5i 

Mdm2 with p53 complex 

X4 

4. 

Mdm2jmRNA 

Mdm2 messenger mRNA 

X3 

5. 

PML 

Unbounded PML protein 

Xh 

6. 

Mdm2.PML 

Mdm2 PML Complex 

xe 

7. 

p53.PML 

p53_PML complex 

X7 


Tables 
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Table 2 List of chemical reaction, kinetic law and their rate constant 


S.No 

Reaction 

N ameoftheprocess 

KineticLaw 

RateConstant 

References 

1 

^'1 1 

X4 - > X4 + X2 

Mdm2 synthesis 

ki{XA) 

4.95 X lO-'^sec-i 

lain]. 

2 

^2 I 

Xi - Xi + X4 

Mdm2jmRNA syn¬ 
thesis 

k2{Xi) 

1.0 X lO-^sec-i 

nig. 

3 

X4 - > (p 

Mdm2jmRNA decay 

ksixi) 

1.0 X 10-Sec-i 

lag. 

4 

A 

X2 ->• <P 

Mdni2 degradation 

k4{x2) 

4.33 X lO-'^sec-l 

lag. 

5 

1 ^5. 

0 -^ Xi 

p53 synthesis 

fcs 

0.078sec“^ 

lag. 

6 

kg 

X3 -)► X2 

Mdm2.p53 degrada¬ 
tion 

keixs) 

8.25 X lO-'^sec-i 

gg. 

7 

Xi + X2 -^ ^3 

Mdm2jp53 synthesis 

kT{xi)(x2) 

11.55 X lO-^sec"^ 

0- 

8 

kg 

X3 ->• Xi + X2 

Mdm2.p5S dissocia¬ 
tion 

ksixs) 

11.55 X 10-®sec-i 

0- 

9 

1 kpM,L 

0 -^ X5 

PML creation 

kpML 

5.0 X lO-^sec-i 

m- 

10 

1 ^'10, 

X2 + X3 - Xq 

PML_Mdm2 complex 
formation 

klo{x2){X5) 

2.0 X 10“^sec“^ 

m- 

11 

fclL 

Xq f x^ 

Decreasing of Mdm2 

fcll(xi) 

3.3 X lO-^sec-^ 

m- 

12 

X5 (p 

PML decay 

kl2{x3) 

5.0 X lO-^sec-i 

g. 

13 

Xi + X5 -^ X7 

p53-PML complex 

formation 

kl3{xi){X3) 

1.0 X 10“®sec“^ 

gg. 

14 

fcl4 

X7 —^ Xi + ^5 

dissociation of 

pbi-PML complex 

kl4{x7) 

1.0 X 10“®sec“^ 

gg. 























Figure Legends 

Figure 1 - PML induced p53 — Mdm2 biochemical regulatory model 

The schematic diagram of stress p53-Mdm2-PML model. 
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Figure 2a - Effect of PML on p53 protein 

Temporal behaviour of p53 due to the switching of PML protein at various concentration level i.e 
0.04,0.01,0.009,0.006,0.003,0.001,0.0005,0.0008. 

Figure 2b - Effect of PML on Mdm2 protein 

Temporal behaviour of Mdm2 due to the switching of PML protein at various concentration level i.e 
0.04,0.01,0.009,0.006,0.003,0.001,0.0005,0.0008. 

Figure 3 - Two dimensional recurrence plot 

Two dimensional recurrence plots between {p53 — Mdm2), {p53 — PML), {Mdm2 — PML), (p53_Mdm2 — 
PML), {Mdm2_mRNA — PML)ioT different values of rate constants k_PML, i.e. 0.0005, 0.001, 0.003, 
0.006, 0.009, 0.01 and 0.03 respectively. 

Figure 4 - Stability Analysis 

Plots for the stability analysis.In the upper panel stability of p53(in magenta) and Mdm2(in blue) is shown 
in terms of variablity in amplitude with respect to time. In lower panel stability of p53(in magenta) and 
Mdm2(in blue) is shown in terms of variablity in time period with respect to time. 
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